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Creation of quantum correlations between two atoms in a dissipative environment from an 
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We have investigated the effect of counter-rotating terms on the dynamics of entanglement and 
quantum discord between two identical atoms interacting with a lossy single mode cavity field for 
a system initially in a vacuum state. The counter-rotating terms are found to lead to steady states in 
the long time limit which can have high quantum discord, but have no entanglement. The effect of 
cavity decay rate on this steady state quantum discord has been also investigated, surprisingly, the 
increase in cavity decay rate is found to both enhance and maximize the steady quantum discord for 
separable states. 
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I. INTRODUCTION 

Rabi model is one of the most complete quantum 
mechanical models used to describe the interaction be- 
tween a bosonic field and a two-level system fT], such 
as a two-level atomic system interacting with an elec- 
tromagnetic field in a cavity in quantum optics ijH], an 
electronic spin coupled to phonon modes of a crystal 
lattice [3], superconducting qubits interacting with a 
nanomechanical resonator [4], and many more exam- 
ples can be found in Ref. [5^. The Rabi Hamiltonian 
reads {fi = 1) 

H=^o-z + coa^a + g{a+ + a-){a + a'^), (1) 

where co and loq are the field and atomic transition fre- 
quencies, respectively, g is the coupling constant, a (a^) 
is the annihilation (creation) operator of the bosonic 
field and <7z, cr+ and a- are the pseudo-spin operators. 
It has been difficult to find an analytic solution to the 
Rabi model. Recently, Braak has shown that its eigen- 
values can be calculated as the roots of a polynomial 
However, a large number of studies on the model in- 
volve numerical perturbative |@] and approximate 
analytical solutions 13] ■ The most widely used ap- 
proximation is the so called rotating wave approxima- 
tion (RWA) which amounts to ignoring the counter- 
rotating terms, a'^a^ and ad-, in the field-atom interac- 
tion [9]. The ignored terms correspond to the emission 
and absorption of virtual photons without energy con- 
servation. The error introduced by RWA depends on 
the magnitudes of the atom-field frequency detuning 
and the atom-field interaction strength. This approxi- 
mation is thought to be valid only for small atom-field 
frequency detuning (i.e., |A| = |a;o — << (^,<^o) ^rid 
weak couplings (i.e., g << co,coo)- Iri general, RWA is 
justified in typical optical setups, because the atom-field 
coupling strength divided by the field frequency is of 
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the order 10"'' — 10~^ and the nearly resonant condition 
can be satisfied fld^. In this limit, Rabi model is known 
as Jaynes-Cummings model which can be integrated ex- 
actly, so it is attractive and has been used successfully 
for more than four decades to explain many physical 
phenornena in quantum optics, such as Rabi oscilla- 
tions [11], squeezing 1IT2I1 , non-classical states (such as 
Schrodinger cat-like states) Q, Fock states col- 
lapse and revival of atomic inversion j l^ and quantum 
as well as classical correlations between atom-atom or 
atom-field systems I1I6I . However, recent developments 
in physical implementation of qubits lead to systems 
with g/co up to one and possibly much higher than 
one which requires a careful consideration of the effect 
of virtual processes (counter-rotating terms) on atom- 
field interaction 10^. 

The quantum correlations are crucial to many quan- 
tum information tasks; for example entanglement is im- 
portant in speeding up quantum algorithms, making 
quantum teleportation [17] and cryptography II8II pos- 
sible. On the other hand, quantum discord (QD) lll9ll 
has a role in the deterministic quantum computation 
with one pure qubit [|20il which was also demonstrated 
experimentally in recent times jzill and in Grover search 
algorithm [22]. Dakic et al. recently showed that QD is 
an optimal resource for remote quantum state prepara- 
tion 12311 . However, the role of QD in speeding up quan- 
tum operations is still controversial. On the other hand, 
all realistic quantum systems are inevitably connected 
to decohering and/or dissipative environments which 
wipes out the quantumness of the systems. The ef- 
fect of decoherence and dissipation on the dynamics of 
quantum discord and entanglement between two qubits 
in contact with an environment have been investigated 
very throughly by many groups in the last decade ex- 
perimentally as well as theoretically fid, '24'- 36]. It 
was shown that the entanglement is very fragile and 
can cease to exist in a finite time, although the coher- 
ence of the single qubit decays exponentially. This ef- 
fect is named as entanglement sudden death (ESD) ll23l 
and also observed experimentally | |25|| . In contrast to 
entanglement, QD is more robust under decoherence 
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and dissipation and it is considered as a measure of 
more-general-than-entanglement type quantum corre- 
lations. Actually, it was shown by Ferraro et al. that 
almost all quantum states of a bipartite system have 
non-zero quantum discord IIstII . It was found that 
QD presents an instantaneous disappearance at some 
time points in non-Markovian regime (2^ and asymp- 
totic decay in Markovian regime even at finite temper- 
atures I2 7[,[2§], w hile entanglement dies in a finite time. 
In RefsjMlssll. QD is found to be unaffected by de- 
coherence or dissipation for certain initial states con- 
trary to the entanglemen t dy namics which suffers sud- 
den death. In Refs. |29l l36ll , quantum correlations, as 
quantified by quantum discord, are found to increase 
under memoryless dissipation for some initially sepa- 
rable states, although the qubits remain unentangled 
for all time. 



The counter-rotating terms are found to be respon- 
sible for several novel quantum mechanical effects, for 
example generation of photons as well as entanglement 
from states with no initial excitation (vacuum) is such 
subject which has received plenty of attention in recent 
years [38-42]. Along these lines, Ficek et al. showed 
that it is possible to create high degree of entangle- 
ment between two qubits in a leaky cavity with the 
help of the virtual processes jS^. An efficient method 
to generate entanglement between two separate ensem- 
bles of molecules trapped inside a superconducting res- 
onator through which they are strongly coupled to a 
microwave field mode has been outlined in Ref. |40|] . 
The possibility of creation of photons as well as atom- 
cavity entangled states in the experimental implemen- 
tation of the Landau-Zener sweeps |41] and the non- 
stationary circuit QED [42] was proposed recently. The 
creation of more-general-than-entanglement type quan- 
tum correlations, such as quantum discord, from ini- 
tially vacuum states is an interesting and significant 
topic and has not received enough attention. To the 
best of our knowledge, there is no reported study con- 
cerning the creation of QD from initially vacuum states. 



In this Letter, we have studied the creation of entan- 
glement and quantum discord between two qubits in- 
teracting with a lossy single mode cavity field and ini- 
tially in a state with zero excitation. We have demon- 
strated that the counter-rotating terms in the Rabi 
Hamiltonian lead to steady states in the long-time limit 
which have no or negligibly small entanglement but can 
have high quantum discord. It is natural to expect that 
the cavity decay decreases the magnitude of the steady 
state quantum discord, but the results we have obtained 
indicate that an increase in cavity decay rate increases 
and also maximizes the magnitude of the steady state 
quantum discord for separable states when the virtual 
processes are taken into account. 



II. THE MODEL 

Here, we consider two noninteracting two-level 
atoms, called A and B, coupled to a single mode cav- 
ity field in a leaky cavity. The dynamics of the atoms 
and the cavity field is determined by the master equa- 
tion for the density operator pg of the atoms plus cavity 
field (for h = l) HHMIH: 

Ps = -i[H,ps] - ^(flVs -2apsa^ +psa^a), (2) 

where k is the damping rate of the cavity mode and 
determines the quality of the cavity with the relation 
Q = lo/k, where co is the field frequency. H in Eq. Q 
is the Hamiltonian of the atom-cavity system includ- 
ing counter-rotating terms and can be obtained by re- 
placing (Tz, cr^ and (T- with (cr^ + a^), [erf + cr^) and 
{cr^ + a^) in Eq. respectively. Here, we have as- 
sumed that the atoms have identical transition frequen- 
cies, ojQ, and coupled to the field mode with the same 
coupling constant, g. 

The interaction terms in the Hamiltonian for two 
atoms interacting with a single-mode cavity field con- 
tain terms such as a'^a''j_ and aaL (for j = A, B) which 
are called coimter-rotating terms. These terms do not 
conserve the total number of excitations in the system: 

the terms determined by aa[_ describe the process in 
which a photon is annihilated in the cavity mode as the 
atom makes a downward transition, while the terms 
a'^a'j^ signify that a photon is created in the cavity mode 
as the atom makes an upward transition. If the counter 

rotating terms, fl^u;'^ and acL for / = A, B, are ignored 
(i.e., when the RWA is made), the Hamiltonian takes 
the form 

+g{ala + a^a^). (3) 

The RWA is valid for small detuning, i.e., |A| = jwo ~ 
oj\ << ojQ,co and for weak couplings, g « CjOq,oj. 

The master equation Q has no analytical solution 
neither for RWA nor for non-RWA Hamiltonians js^ 
H^, but can be solved numerically for different initial 
states. In the present work, we will restrict ourselves 
to the pure product state with zero initial excitation, 
i.e.,ps(0) = \gA,gB,^) {SA,gBA- Here \gA,gB) denotes 
that the atoms, A and B, are in their ground states and 
|0) signifies that there is no photon inside the cavity. 



III. CORRELATION MEASURES: ENTANGLEMENT 
AND QUANTUM DISCORD 

The reduced density matrix of the atoms A and B 
can be calculated by tracing the density operator, p^, 
over the cavity degrees of freedom. For the initial 
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state considered in the present work, the reduced den- 
sity matrix of the atoms in the two-qubit standard ba- 
sis |1) = |e^,eB>,|2) = |eA/^B>,|3) = \gA,eB) and 
|4) = IgAfSs) has the X structure: 



Pab 



Pn pu 

,022 P23 

P32 P33 

,041 



(4) 



Based on symmetry considerations, it can be shown that 
the dynamics in Eq. ^ preserves the X-form of the den- 
sity matrix O . By using the density matrix the 
quantum correlations between the atoms as measured 
by entanglement and quantum discord can be calcu- 
lated. We adopt Wootters' concurrence as entan- 
glement measure which is a normalized measure of en- 
tanglement and gives for separable states and 1 for 
maximally entangled (Bell) states. For the density ma- 
trix l(4|, the concurrence fimction reads 

C(f) = 2max{0,Ci(0,C2(0}, 

Cl(0 = ||023| - VPnPU' ^2(0 = ||0i4| - ^pnP33- 

(5) 

On the other hand, quantum discord captures non- 
classical correlations between two two-level systems 
that are more general than entanglement IT9I1 . It is de- 
fined as 



D{t) = m-](t), 



(6) 



where /(f) = S{p a) + S(|Db) — '^'{Pab) is the total cor- 
relations between the atoms; S(^p) = —Tr{plogp) is the 
von Neumann entropy and (jOg) is the reduced den- 
sity matrix obtained by tracing p^^g over the subsys- 
tem B (A). The other quantity }{t) is the measure of 
classical correlations between the atoms defined as the 
maximum information one can get about the atom A 
(or B) by performing a set of von Neumann type mea- 
surements on the atom B (or A) [27]. Obviously, }{t) 
and also QD are not symmetric quantities; i.e., they de- 
pend on which the measurement is performed. Here 
we shall consider a set of positive-operator-valued mea- 
surements performed on the subsystem B | 46] . Re- 
cently, analytical expressions fo r Q D of X-state density 
matrix have been reported ll46l - l48ll . Here, we will use 
the results given in Ref. m which are, in fact, equiva- 
lent to that of Ali et al. given in Ref. jiTll . According to 
the results in Ref. fiill , quantum discord is given as 

D(f) =min{Qi,Q2}, (7) 

where Qj = M{pn + P33) + Ef=i A^g logj Aj^g + Py, 
with A^g being the eigenvalues of pABr = 

M(t), P2 = -EtlP/7log2Pn - + 

T = (1+ V(l-2(p33 + P44))^ + 4(|p23| + |pi4|)^)/2 

and M(a) = — alogjft: — (1 — a)log2(l — a) is the 
binary Shannon entropy function. 



For pure states, entanglement of formation and quan- 
tum discord are found to be equivalent, while for mixed 
states such an identification is much more difficult to 
make. These two measures can disagree on the quan- 
tum correlations of a mixed state, for example, it was 
shown that some unentangled mixed states can carry 
non-zero quantum discord |i33]. One should note that 
although QD measures non-classical correlations that 
are more general than entanglement, it is not a faithful 
non-classical correlation measure since it does not van- 
ish only for a state which is strictly classical correlated 
(see Ref. iH). 



IV. RESULTS 

In the following, we will investigate the creation of 
entanglement and quantum discord between two atoms 
interacting with a lossy single mode cavity field and ini- 
tially in a state with zero excitation. To do this, we will 
solve the master equation l|2]l for the non-RWA Hamil- 
tonian and we will use Eqs. (|5) and ^ to calculate con- 
currence and QD, respectively. In our calculations, we 
will fix the detuning K = loq — lo = Q.OIlo. 

A brief outline of the numerical procedure to solve 
Eq. @ should be given. Due to excitation number be- 
ing a conserved quantity in Jaynes-Cummings model, 
the Hamiltonian is block-diagonal in the Hilbert space 
of ® (g) F~ (where and F~ indicate the qubit 
and 00-dimensional Fock space for the field, respec- 
tively), and diagonalizrng blocks corresponding to to- 
tal excitations of the system gives the analytic solu- 
tion of the problem. In the full-model, excitation con- 
servation is no longer valid and one should diagonal- 
ize an infinite dimensional Hamiltonian. In practice, 
for the considered initial state we have done a conver- 
gence study of the Fock space dimension of the cav- 
ity field filMlillli. h\ this work, we have consid- 
ered basis vectors of type \iA/jB/ n) where i,] = e,g and 
n = 0, 1,2, N. The converge criterion was consid- 
ered as the absolute value of the system density matrix 
elements was smaller than 10"^*^ which was found to 
be satisfied approximately for N = 38 for the largest 
atom-field coupling constant and the smallest cavity de- 
cay considered in the present work; for the results re- 
ported in the remainder of the text, we have taken into 
account all the basis vectors |fA/7B/W) where = e,g 
and n = 0, 1,2, ...,50. Also, for all the considered cases 
the basic properties of the atom-atom density matrix, 
such as positivity, hermiticity and trace preservation, 
have been checked during the computational process. 

First, we investigate the effect of atom-field cou- 
pling strength on the time evolution of concurrence and 
quantum discord. Fig. 1 displays the time-dependent 
QD and concurrence for the initial state ps{0) = 
\gA,gB,0) {gA,gB,0\ with A = (a;o - a;) = O.Olo; and 
K = 0.2a;. Fig. 1(a) is for the so called weak-coupling 
regime where g = 1.0 x 10~*a;, g = 7.5 x lO^^o^, 



4 



(a) 



0.008 



0.002 











1-1.5x10- 






g 3 5x10" 


s 




.-g-2.5xl0- 













5.0x10' 1.0x10" 1.5x10* 2.0x10" 2.5x10" 3.0x10" 




100 120 



FIG. 1. Effect of atom-field coupling strength, g, on the dy- 
namics of quantum discord and concurrence (inset in (b)) ver- 
sus cvt for ps{0) = \gA,gB,0) {gA,gB,0\ initial state, cvq = 
LOlw, K = 0.2lo and for two-qubit Rabi model. Here (a) 
describes the weak coupling regime and includes plots for 
^ = 1.0 X lO^^'a; (red, dashed), g = 7.5x IQ-'^iv ( green, solid). 



^ = 3.5 X lO^^a; (blue, dotted) and g = 2.5 x lO'^cv (black, 
dot-dashed), while (b) indicates strong coupling regime and 
includes plots for g = cv (red, dashed), g = 0.75a; (green, 
sohd), g = 0.35a; (blue, dotted) and g = 0.25ai (black, dot- 
dashed). Note that for weak coupling no entanglement is in- 
duced, so they are not plotted here. 



^ = 3.5 X lO^^o; and ^ = 2.5 x lO^^o;, while Fig. 1(b) is 
for the strong coupling regime with g = g = 0.75a;, 
g = 0.35a; and g = 0.25a;. It is expected that the 
difference between RWA and non-RWA evolutions will 
be more pronounced for the strong coupling parame- 
ters. Since RWA dynamics conserve the total excitation 
number and the considered initial state has zero exci- 
tation to start with, both QD and concurrence remain 
zero at all times. On the other hand, counter-rotating 

terms (fl^cr^ and aaL for j = A, B) will produce virtual 
excitations which create quantum correlations between 
the atoms as can be seen from Figs. 1(a) and 1(b). Even 
in the weak coupling regime, one can see the effect of 
counter-rotating terms as displayed in Fig. 1(a); non- 
RWA d5mamics lead to steady-state non-zero quantum 
discord (D = 0.0025), albeit very small. The magni- 
tude of ^ in this regime seems to control the speed of 
approaching the steady-state QD value rather than the 
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FIG. 2. Time evolution of the mean excitation number in 
the whole system, {Nj)^ (red-dashed and black-solid), and 
in the atomic system, (Nyi)j (blue-dotted), versus cvt for 
Ps{0) = \gA,gB,0) {gA,gB,0\ initial state, wq = l.Olw, g = cv 
and K = (red-dashed) and k = 0.2a; (black-solid and blue- 
dotted). Here, (Nj), = {a'^a + a^a^ + alo-^)^ and (N^), = 
(^A^A + ^B^B^^ ^here (. . = Tr {ps{.. .}). 



magnitude of the asymptotic QD. The dynamics of QD 
and concurrence in the strong-coupling regime under 
non-RWA dynamics is richer than the one imder the 
weak coupling; both concurrence and QD are found 
to be non-zero as displayed in Fig. 1(b). The induced 
entanglement goes through a series of sudden deaths 
and births and finally dies out with an overall lifetime 
inversely proportional to the coupling constant g. In 
contrast, quantum discord always approaches a non- 
zero asymptotic value which has no monotonic relation 
with g. The observed behavior of QD suggests that 
not all quantum correlations created by virtual excita- 
tions are lost in the dissipative d5mamics of the atom- 
cavity system. Note that QD approaches its steady 
value much faster for strong coupling than for weak 
coupling. In Fig. 2, we display the average excitation 
number {Nt)^ = (^a^a + a^a^ + cr^cr't)^ in the system 
as a function of dimensionless time, cot, at g = to for 
K = and K = 0.2a;, along with average excitation 
number, (Ny^)f = {(rfcr^ +cr^cr^)f, in the atomic sub- 
system for g = LO and k = 0.2a;. Due to the collective 
radiation inhibition effects, the so called "virtual pho- 
tons" produced by the counter-rotating terms remain in 
the cavity despite the strong dissipation |50] (see solid 
line in Fig. 2). Fig. 2 also shows that the mean exci- 
tation number of the atomic subsystem for leaky cav- 
ity is also frozen in the long time limit which is nearly 
1. This leads to the appearance of steady states which 
can have high QD. Moreover, as can be seen from this 
figure, {Nt)i is non-zero for both k = 0.2a; and the 
non-dissipative cavity (k = 0). Non-zero (Nj)^ is some- 
times claimed to be due to the dissipation inhibiting 
the destruction of virtually created photons in the cav- 
ity issllsill which seems to be not the case because even 
for K = 0, {NT)f has a high value. 

Now, we consider the effect of cavity decay on the 
dynamics of concurrence and QD for the same ini- 
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FIG. 3. Effect of cavity decay on the dynamics of quan- 
tum discord and concurrence (inset in (b)) versus ivt for 
Ps{0) = \gA,gB,Q) {gA,gB,0\ initial state, coq = l-Olw, k = 
0.08w (black, dot-dashed), k = 0.2a; (blue, dotted), k = 
2a; (green, solid), k = 20a' (red, dashed) and without the 
RWA. Here (a) describes the weak coupling regime and in- 
cludes plots for ^ = 3.5 X lO^^o;, while (b) indicates strong 
coupling regime and includes plots for g = 0.35a;. Note that 
for g = 3.5 X lO^^o; and all considered cavity decay rates as 
well as for g = 0.35a; and k = 20a;, no entanglement is in- 
duced, thus they are not plotted here. The inset in (a) shows 
QD versus cot for small QD region. 



tial state {\gA,SB^O)) and detuning (A = O.Olo;) in 
weak ( g = 3.5 X 10~^a;) and strong {g = 0.35a;) cou- 
pling regimes. We display the QD and concurrence at 
K = 0.08a;, k = 0.2a;, k = 2a> and k = 20a; for weak 
and strong coupling regimes in Figs. 3(a) and 3(b), re- 
spectively. Increasing the cavity decay rate is found to 
decrease the maximum of induced entanglement and 
to hasten its death at strong coupling regime, as ex- 
pected (see the inset in Fig. 3(b)). On the other hand, 
the decay rate dependence of quantum discord is sur- 
prising for both the weak and strong coupling param- 
eters as can be seen from Figs. 3(a) and 3(b); the in- 
creasing the cavity decay rate is found to enhance the 
steady state QD. For high values of k (k = 2a; and 
K = 20a;) QD approaches relatively high steady-state 
value independently of whether g is large or small, 
while for lesser decay rates (k = 0.08a; and k = 0.2a;) 



the steady-state QD is quite different for the weak and 
strong coupling cases; the steady state QD is promi- 
nently high for g = 0.35a; compared to the one for 
^ = 3.5 X lO^^a;. One should note that the steady 
state is reached in much longer times for the weak- 
coupling case compared to that of the strong-coupling. 
These findings can be explained based on the results in 
Ref. iSQ] where the author studied the effects of counter 
rotating terms on the expectation value of the popu- 
lation inversion operator (energy shift) for N identical 
two level atoms interacting with a single mode cav- 
ity field in a leaky cavity (dissipative Dicke Model). 
It was found that in very bad-quality cavities, signifi- 
cant steady-state energy shifts can be obtained which 
increases with the increase in k, while these shifts do 
not depend on the size of the atom-field coupling con- 
stant. On the other hand, it was also demonstrated that 
these energy shifts are small and mainly depend on g 
for high-quality cavities. From the above results, the 
large cavity damping is found to increase the contribu- 
tion of counter-rotating terms and provides us a way 
to see the effects of virtual processes at low atom-field 
coupling constants. The role of k in the creation of long 
time non-zero QD can be better understood by look- 
ing at the atom-atom density matrix; for example, for 
A = O.Ola;, g = 0.35a; and k = 0.08a;, the steady state 
density matrix is equal to 



Pab{k = 0.08a;) 



0.014 

0.073 

0.073 

0.074 



0.074 ■ 
0.073 
0.073 

0.84 



with steady discord D = 0.1, while for k = 20a;, 



Pab{k = 20a;) 



0.32 



\ 3.1 X 10-5 



3.1 X 10-5 \ 

0.17 0.17 

0.17 0.17 

0.34 J 



with steady discord D = 0.33. Comparing these 
two density matrices one can note that the popula- 
tions pn = \eA,eB) {eA,eB\, pii = \eA,gB) {eA,gB\ 
and ,033 = \gA,eB) {g A, sbI and the coherence 
P23 = \^A/gB) {gAf^Bl increase with the increase 
of K, while the population ^44 = \gA,gB) {gA,gB\ and 
the coherence pu = |eA/^B) {gA>gB\ decrease. Such a 
change can enhance the steady state quantum discord. 
One can conclude that the increase in the cavity decay 
can decrease the chance of atoms to interact with cavity 
field and the atomic system is frozen in a state with 
high quantum correlations before it has a chance to 
decay to the groimd state Issll . It is interesting to note 
that in the very bad quality cavity case {k/lo > 20), the 
steady atomic density matrix is saturated and can be 
nearly written as an equal weighted sum, pAg{K/co > 
20) = 0.333 {\eA,eB){eA, eel + \g A, gB) {gA, gB\) + 
0.334 (T+|, where = (|e^,gB> + \gA,eB)) 
is the Bell state, and the ^-dependence of this atomic 
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steady state becomes negligible as shown in Fig. 3. 
Although this atomic state contains no entanglement, 
it has high QD which is nearly D = 0.333. Recently, 
it was shown that maximal QD reachable by two 
separable qubits is 1/3 L52i] . As a consequence, it 
seems that the interplay between losses and virtual 
processes does not only lead to high steady QD, but 
also maximizes it for separable states. A similar result 
has been also obtained in different systems (see, for 
example Ref. [46]). One should note that a negligibly 
small amount of steady entanglement (C ~ 2 x 10^'') 
exists in the state only for k = O.OSci; that can be noted 
from the above reduced density matrix. 

V. CONCLUSION 

We have investigated the effect of counter-rotating in- 
teraction terms on the dynamics of entanglement and 
quantum discord between two qubits interacting with a 
single mode cavity field in a leaky cavity for a system 
with zero initial excitation. We have shown that virtual 
processes can lead to long-lived constant quantum dis- 
cord even at weak atom-field couplings {g/to k 10~^) 
and can create non-zero entanglement at strong cou- 
pling regime which goes through sudden death and 
birth processes before permanently vanishing (for k > 
0.08a;). 

A counter-intuitive finding of the present work is 
the cavity decay rate dependence of entanglement and 



quantum discord; while an increase in decay rate leads 
to a shorter lifetime for entanglement, it increases the 
steady state value of quantum discord for both strong 
and weak couplings. Moreover, for high quality cav- 
ities (k = 0.08a; and k = 0.2a;), the steady state QD 
is found to be highly interaction strength dependent, 
while for low-quality cavities (k = 2a; and k = 20a;), 
the dependence of steady state QD on the interaction 
strength is found to be negligible. In fact, the competi- 
tion between coimter-rotating terms and cavity decay is 
found to give not only a relatively high value of steady 
QD, but indeed maximizes it for separable states in the 
case of very bad quality cavity (k/co > 20). 

It is worth mentioning here that the master equa- 
tion (|2) has been used by a large number of groups 
in recent years in order to study the cavity decay in 
Rabi model [3^111,141,113], but its validity has not been 
proved yet [3§\. Nevertheless, if the master equation 
is indeed applicable, the results reported in this Letter 
might be relevant for the strong coupling experimen- 
tal work which is madepossible with the recent cavity 
QED circuit proposals 10-111 • 
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